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Abstract 


The QED contribution to the anomalous magnetic moments of electron and muon are known 
very precisely up to the order a'^. However, the knowledge of term will also be required when 
the precision of measurement improves further. This paper reports the first systematic attempt 
to evaluate the term. Feynman diagrams contributing to this term can be classified into six 
gauge-invariant sets which can be subdivided further into 32 gauge-invariant subsets. Thus far 
we have numerically evaluated all integrals of 17 gauge-invariant subsets which contain light-by- 
light-scattering subdiagrams and/or vacuum-polarization subdiagrams. They cover most of leading 
terms of muon g — 2 and lead to a preliminary result 663 (20) (a/vr)^, which is 8.5 times more 
precise than the old estimate. 
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I. INTRODUCTION 


The deviation of the electron g value from 2 predicted by Dirac’s theory was first con- 
hrmed by an experiment on atomic spectrum Schwinger showed that this deviation can 
be explained as the effect of radiative correction by the relativistic renormalized QED which 
he had developed [3]. Together with the discovery of Lamb shift in the spectrum of hydrogen 
atom ^ , it provided convincing experimental evidence that (until then discredited) QED is 
capable of predicting the effect of electromagnetic interaction precisely, provided that it is 
renormalized. 


A. Measurement of electron g-2 


By 1970’s the precision of measurement of electron g-2 was improved by four more orders 
of magnitude by means of spin precession of the electron moving in a constant uniform 
magnetic held j^. The value of the electron g-2 was improved further by three additional 
orders of magnitude in a Penning trap experiment by Dehmelt’s group at the University of 
Washington. Their published results are p 

ae- = 1 159 652 188.4 (4.3) x lO’^^ [g g 

ae+ = 1 159 652 187.9 (4.3) x lO"^^ [3.8 ppb], (1) 


where the numerals 4.3 in parentheses represent the combined statistical and systematic 
uncertainties in the last digits of the measured value. 1 ppb = 10“®. 

The precision of measurement has thus been improved by seven orders of magnitude over 
40 years. This enormous improvement in measurement was matched by the improvement of 
theory of radiative correction to the electron g-2 from the order a to the order a^, leading 
to the most stringent test of the validity of QED. 

The uncertainty of the experiment (CD was dominated by the cavity shift due to the inter¬ 
action of the electron with the hyperboloid cavity, which has a veiw complicated resonance 
structure. Several efforts were made to reduced this uncertainty p, Qj. One of them is to 
replace the hyperboloid cavity W a cylindrical cavity, which allows analytic computation of 
the structure of the resonance • Gabrielse’s new measurement of the electron g-2 is based 
on this analysis. Recently a preliminary result of this measurement was reported, which is 
7.5 times more precise than cdS 
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B. Theory of electron g-2 to order 


( 2 ) 


(3) 


The QED contribution to the electron g-2 can be written as 

tteiQED) = Ai + A2{me/m^) + A2{me/mr) + A2,{nie/m^,, me/m^) 
and Aj, i = 1, 2, 3, can be expanded as 

The first four coefficients of Ai are 

= 0.5, 

Af^ = -0.328 478 965 ..., 

Af^ = 1.181 241 456..., 

= -1.728 3 (35). 

A^^'^ and A^f’ are known analytically was obtained by both numerical lid 

and analytic integrations is obtained thus far by numerical integration only jl^ . 

Its uncertainty has been reduced by 10 compared with the old one j^. Although it has 
been evaluated by one method only, it has been subjected to an extensive cross-checking 
among diagrams of 8 th-order and also with 6 th-, 4th-, and 2nd-order diagrams. 

A 2 terms are small: 

A^f^{me/rn^){a/TiY = 2.804 x 10“^^, 

A^f^{me/rnr){cilTiY = 0.010 x 10 “^^, 

A 2 \f^e/‘m^){a/'K)^ = —0.924 x 10“^^, 

Af\me/rnr){a/'K)^ = —0.825 x 10“^® (5) 

ler (~ 2.4 x 10“^^). The non-QED contribution of 

1171, Ui^ 


The contribution of A 3 term is even sma 
the Standard Model are also known ll 


ae(hadron) = 1.671 (19) x 10 
ae(weak) = 0.030 (1) x 10“^^. 


( 6 ) 


To compare the theory with the measured value of Og one needs a value of a obtained by 
some non-QED measurement. The best available a at present is 3,0] 


a-\h/Mc,) = 137.036 000 1 (11) [7.4 pph] 


(7) 
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This leads to 


ae{h/Mc,) = 1 159 652 175.86 (0.10)(0.26)(8.48) x 10"^^ 
ae{exp) — ae{h/McJ = 12.4 (9.5) x 10“^^, 


( 8 ) 


where 0.10 is the remaining uncertainty of the term, 0.26 is based on an educated guess 
= 0(3.8)) made by Mohr and Taylor and 8.48 is the uncertainty in the measure¬ 
ment ©. The error 8.48 is still large but is within a factor 2 of the error of the Seattle 
measurement (P). 


C. Measurement of muon g-2 


The last and best of three measurements of the muon g-2 at CERN had an uncertainty 


of 7 ppm j2]J. After years of hard work the muon g-2 measurement at the Brookhaven 
National Laboratory has come close to the design goal (0.35 ppm): 221 


a^-(exp) = 11 659 214 (8) (3) x 10 


-10 


(0.7 ppm). 


(9) 


The world average of a^(exp) obtained by combining this and earlier measurements 


24, 3 is 


a^(exp) = 11 659 208 (6) x 10 


-10 


(0.5 ppm). 


( 10 ) 


D. Hadronic and electroweak contributions to muon g-2 


Currently, the prediction o 
of the hadronic contribution 


le Standard Moc 


26 


27 


281 


3G, 


el reflects the difficulty in the treatment 


31|. The lowest-order hadronic vacuum- 
polarization effect on has thus far been determined from three sources, 

(i) e’''e“ annihilation cross section, 

(ii) hadronic r decays, 

(iii) e’''e“ —> 7 - 1 - hadrons [radiative return]. 

The process (i) has been analyzed by many groups over years. Some recent results are 




33|,|3i 


a^(had.LO) = 6934 (53)^.^ (36)™d x 10-i\ 
a^(had.LO) = 6924 (59)^.^ (24),,^ x 10-i\ 
a^(had.LO) = 6944 (48)exp (10)™d x IQ-'b 


( 11 ) 
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Recent estimate of hadronic light-by-light scattering contribution is ^ 

a^(had.NL) = 136 (25) x 10“^^. 

The electroweak interaction effect is known to two-loop order: l^| 

a^{weak) = 152 (1) x 10“^^, 
a^{weak) = 154 (1) (2) x 10“^^, 


( 12 ) 


(13) 


where ( 1 ) and ( 2 ) in the second line are estimates of remaining theoretical uncertainty and 
Higgs mass uncertainty, respectively. 

E. QED contribution to muon g-2 to order 


The QED contribution to the muon g-2 can be written as 


a^iQED) = Ai + H 2 (m^/me) -h H 2 (m^/me) + H 3 (m^/me, m^/me). 


(14) 


Renormalizability of QED guarantees that the functions Ai, A 2 , and A^ can be expanded 
in power series in o/tt with finite calculable coefficients: 

2 /-N /q;\3 


A, = H 


f. (^) , .y (^) + 4. + ,,.. , . , 3. 


( 16 ) 


Ai in (HI is identical with Ai of m and has been evaluated to the eighth (i.e., order. As 
for A 2 and A 3 , it is easy to see that A^^'^ = A^'^ = A^^'^ = 0 since they have no corresponding 
Feynman diagram. A^\ A'^\ and terms have been evaluated accurately by power series 
expansion in me/m^ or me/m-r. Thus far A^^^ and A® are known mostly by numerical 
integration. 

The QED contribution a^(QED), even though it is the predominant term of a^, has 
received little attention for many years because of its small error bars. The theoretical 
uncertainty came predominantly from the term whose contribution to is about 3.3 
ppm. Recently we have completed a new evaluation of the term in which all contributing 
terms have been evaluated by two or more independent calculations, uncovering an error in 
some diagrams Q and eliminating the term as a possible source of theoretical uncertainty 
This causes the QED contribution shifted to 


a^(QED) = 116 584 719.43 (0.02)(1.15)(0.85) x 10 


-11 


( 16 ) 
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where 0.03 is the remaining uncertainty of the term, an improvement of factor 40 over the 
previous result Q. The uncertain., 1.15 con.ea Icon, a cade ertin.ale of fhe conrtibufion 
of the term, which now stands out as the largest source of uncertainty in a^(QED). The 
error 0.85 comes from the uncertainty in the value of a given in ( 0 ). 

The prediction of the Standard Model, including the hadronic vacuum-polarization, 
hadronic light-by-light contributions, and the electroweak effect, is 

af,{SM) = 116 591 870.7 (76.2) x (17) 

where the uncertainty in theory is mostly due to the hadronic vacuum-polarization term. 


F. Tenth-order term: Why is it needed ? 

A very important byproduct of the study of is that it gives the best value of the 
hne structure constant a available at present. If we use the new experiment by the Harvard 
group, the precision of a is almost an order of magnitude better than any other measurement 
of a. 

Furthermore, the uncertainty of this measurement is only a factor 2 larger than that of 
theory, which is mostly from the term, since the term is known with small error. Thus, 
when the measurement of is improved further, reduction of the uncertainty of term 
will become crucial in order to obtain a better a(ae). 

For the muon the old estimate of was 930 (170), which contributes only 

0.054 ppm to a^, well within the current experimental uncertainty. Thus improving 
A^ 2 ^\m^/me) is not urgent. However, it will become an important source of error in the 
next generation of experiment. This is why it is desirable to obtain a better value of 
A^ 2 ^\m^/m(.)- The preliminary value of was reported in |^ . 


II. CLASSIFICATION OF TENTH-ORDER DIAGRAMS 


Thus far only a small portion o 


have been evaluated analytically 


tenth-order diagrams contributing to the muon g-2 


4l| |. or numerically Rough estimates based on the 


renormalizati on g roup and other considerations have been made in order to identify leading 


terms 


4al42 


m 


3 - 


6 



















Of course an enormous amount of work is required to go beyond this and evaluate a® 
terms completely. Fortunately, for the muon g-2 the leading contribution comes from those 
Feynman diagrams which contain ln(m^/me) terms whose sources can be readily identihed 
as light-by-light-scattering subdiagrams and vacuum-polarization insertions. Thus relatively 
modest amount of work will enable us to improve the value of over the previous 

crude estimate. On the other hand, the electron g-2, in particular term, is much harder 
to evaluate. Besides its gigantic size none of 12672 diagrams is dominant so that every term 
must be evaluated accurately. 

For both Og and the hrst step is to count and classify Feynman diagrams contributing to 
the term. The contribution to the mass-independent term may be classihed into six 
gauge-invariant sets, which further subdivided into 32 gauge-invariant subsets, depending on 
the nature of subdiagrams (of the vacuum-polarization {v-p) type or light-by-light-scattering 
(/-/) type). Classihcation for A^ 2 ^\m^/me) follows readily from that of With help 

of the Feynman Diagram auto-generator of GRACE system 
diagrams belonging to each set. 


47j | we count the number of 


A. Notation 

Vacuum-polarization functions 11 needed for the evaluation of the contribution of the 
lepton g — 2 may be classihed as follows: In the following denotes a fla, containing 11^ 
on its internal photon line. 

112 , which consists of one closed lepton loop of second-order. 

114 , which consists of three proper lepton loops of fourth order. 

114 ( 2 ); which consists of three diagrams of type 114 in which one 112 is inserted in the internal 
photon line. 

fig, which consists of 15 proper sixth-order lepton loops which do not contain v-p loop. 

116(2), which consists of 30 diagrams of type fig in which one 112 is inserted in the internal 
photon line. 

114 ( 4 ), which consists of 9 diagrams of type 114 whose internal photon line contains a 114 . 
114 ( 2 , 2 ), which consists of 3 diagrams of type 114 whose internal photon line contains two 112’s. 
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rig, which consists of 105 proper eighth-order lepton loops. 


Light-by-light scattering type functions A needed for the evaluation of the contribution 
of the lepton g — 2 may be classified as follows: 

A4, which consists of six proper fourth-order lepton loop. 

( 2 ) 

A 4 , which consists of 60 diagrams in which lepton lines and vertices of A 4 are modified by 
2 nd-order radiative corrections. 

A^'^^ which consists of 105 diagrams in which lepton lines and vertices of A 4 are modified by 
4th-order radiative corrections. 

( 2 ) 

A 4 j^ 2 )) which consists of 60 diagrams in which a v-p diagram 112 is inserted in the photon line 
of Af. 

Ag, which consists of 120 proper lepton loops to which 6 photon lines are attached. 

Finally we need magnetic moment contributions of various orders. 

M 2 , which consists of a lepton vertex diagram of second-order. 

Af 4 , which consists of six proper lepton vertices of fourth-order. 

Mq, which consists of 50 proper lepton vertices of sixth-order. 

Mqll, which consists of 6 sixth-order vertex diagrams containing an external light-by-light 
loop A 4 . 

Mg, which consists of 518 proper lepton vertices of eighth-order which have no closed lepton 
loops. 

f2) 

MgLLb, obtained by replacing A 4 by A 4 in Mqll, and consists of 60 eighth-order vertex 
diagrams. 

MgLLc, obtained by attaching a virtual photon line to the open muon line in all possible 
ways, and consists of 48 eighth-order vertex diagrams. 

MgLLd, obtained by inserting a A 4 internally in the fourth-order M 4 in all possible ways. 




1 (a) 







1(e) 1(f) 





I(J) 


FIG. 1: Set I 


B. Set I 


Diagrams of Set I are built from the magnetic moment contribution M 2 of the second- 
order proper vertex. It consists of 208 Feynman diagrams, which can be classihed further 
into ten gauge-invariant subsets as indicated in Fig. ^ Let us mention only contributions 
to and ^ 2 ^°^ in this and subsequent sections although many subsets contribute also to 
This is because the contribution of ^ 3 ^°^ is negligible compared with others. 

Subset 1(a). Diagrams obtained by inserting four 112 ’s in ^ 2 - One Feynman diagram belong¬ 
ing to this subset contributes to The contribution to comes from 15 diagrams 

(2^ - 1 = 15). 

Subset 1(b). Diagrams obtained by inserting two n 2 ’s and one 114 in M 2 . Nine Feynman 
diagrams of this subset contribute to . The contribution to comes from 63 diagrams 
( 9(23 _ _ g3)_ 

Subset 1(c). Diagrams containing two n 4 ’s in M 2 . There are nine Feynman diagrams that 
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contribute to The contribution to comes from 27 diagrams (9(2^ — 1) = 27). 

Subset 1(d). Diagrams obtained by insertion of one 112 and one 114 ( 2 ) in M 2 . Six Feynman 
diagrams contribute to The contribution to comes from 42 diagrams (6(2^ —1) = 

42). 

Subset 1(e). Diagrams obtained by insertion of one 112 and one Ilg in M 2 . Thirty Feynman 
diagrams contribute to The contribution to comes from 90 diagrams (30(2^ —1) = 

90). 

Subset 1(f). Diagrams obtained by insertion of 114 ( 2 , 2 ) in M 2 . The number of diagrams 
contributing to ^ 4 ^°^ is 3. The contribution to A^ 2 ^^ comes from 21 diagrams (3(2^ —1) = 21). 

Subset 1(g). Diagrams obtained by insertion of 114 ( 4 ) in M 2 . The number of diagrams 
contributing to is 9. The contribution to A^ 2 ^^ comes from 27 diagrams (9(2^ —1) = 27). 

Subset 1(h). Diagrams obtained by insertion of 116 ( 2 ) in M 2 . The number of diagrams 
contributing to is 30. The contribution to A^^^ comes from 90 diagrams (30(2^ — 1) = 
90). 

Subset I(i). Diagrams obtained by insertion of Ilg in M 2 . The number of diagrams con¬ 
tributing to is 105. The contribution to A^^'^ also comes from 105 diagrams. 

Subset I(j). Diagrams obtained by insertion of eighth-order photon propagators, which con¬ 
sist of two A 4 ’s with three photons contracted, in M 2 . The number of diagrams contributing 
to is 6 . The contribution to comes from 18 diagrams (6(2^ — 1) = 18). 

The total number of diagrams of Set I contributing to is 208. The number of 

diagrams of Set I contributing to A^ 2 ^'^ is 498. 


C. Set II 

Diagrams of Set II are built upon six proper fourth-order vertices M 4 by insertion of 
various closed electron loops. It consists of 600 Feynman diagrams, which can be classihed 
further into six gauge-invariant subsets as indicated in Fig. |21 

Subset 11(a). Diagrams obtained by inserting three 112 ’s in M^. The number of diagrams 
contributing to is 24. For the sake of programming convenience this set is subdivided 
into Il(ai), in which all three 112 ’s are inserted in the same photon line, and 11 ( 02 ) in which 
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11(a) 


11(b) 


11(c) 



11(d) 11(e) 11(f) 


FIG. 2; Set II 

112 are on both photon lines. The total contribution to conies from 168 diagrams 

(24(2^ - 1) = 168). 

Subset 11(b). Diagrams obtained by inserting one 112 and one 114 in M 4 . The number of 
diagrams contributing to is 108. The contribution to comes from 324 diagrams 
(108(22 - 1) = 324). 

Subset 11(e). Diagrams obtained by insertion of 114 ( 2 ) in M 4 . The number of diagrams 
contributing to is 36. The contribution to ^ 2 ^'^^ comes from 108 diagrams (36(22 _ _ 

108). 

Subset 11(d). Diagrams obtained by insertion of fig in M 4 . The number of diagrams con¬ 
tributing to is 180. The contribution to comes also from 180 diagrams. 

Subset 11(e). Diagrams obtained by insertion of internal light-by-light diagram A 4 in M 4 . 
The number of diagrams contributing to is 180. The contribution to A^ 2 ^'^ comes also 
from 180 diagrams. 

Subset 11(f). Diagrams obtained by insertion of internal light-by-light diagram A 4 and 
additional 112 in dff 4 . The number of diagrams contributing to is 72. The contribution 
to ^ 2 ^°^ comes from 216 diagrams (72(22 _ _ 2 I 6 ). 

The total number of diagrams of Set 11 contributing to is 600. The number of 

diagrams of Set 11 contributing to is 1176. 
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Ill(a) 


Ill(b) 


III(c) 


FIG. 3; Set III 



FIG. 4: Set IV 


D. Set III 

Diagrams of Set III are built from 50 proper sixth-order vertices Mg by insertion of various 
closed electron loops. This set consists of 1140 Feynman diagrams which can be classihed 
further into three gauge-invariant subsets as indicated in Fig. IHl 

Subset III (a). Diagrams obtained by inserting two 112 ’s in Mq. The number of diagrams 
contributing to is 300. The contribution to comes from 900 diagrams (300(2^ — 
1) = 900). 

Subset Ill(b). Diagrams obtained by inserting 114 in Mg. The number of diagrams contribut¬ 
ing to ^ 4 ^°^ is 450. The contribution to comes also from 450 diagrams. 

Subset III(c). Diagrams obtained by insertion of an internal light-by-light diagram A 4 in 
Mg. The number of diagrams contributing to is 390. The contribution to comes 
also from 390 diagrams. 

The total number of diagrams of Set III contributing to is 1140. The number of 
diagrams of Set II contributing to is 1740. 
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FIG. 5: Set V 


E. Set IV 

Set IV is built from 518 proper eighth-order vertices Mg by inserting a closed lepton loop 
of second order 112 . It has only one subset. Total numbers of diagrams of Set IV contributing 
to and ^ 2 ^°^ are both 2072. A representative diagram is shown in Fig. EJ 


F. Set V 

Set V consists of proper tenth-order vertices. It has only one subset consisting of 6354 
diagrams and contributes only to A representative diagram is shown in Fig. El 

G. Set VI 

This set consists of vertex diagrams of various orders which contain at least one l-l 
subdiagram. Most diagrams of Set VI are built starting from the sixth-order diagram Mqll 
which contains an external light-by-light-scattering subdiagram A 4 . An exception is one 
subset that contains a subdiagram Ag. Note also that diagrams already contained in the 
sets I, II, and III are excluded. The set VI consists of 2298 Feynman diagrams which can 
be subdivided into eleven gauge-invariant subsets as indicated in Fig. El 

Subset VI(a). Diagrams obtained by inserting two n 2 ’s in the internal photon lines of Mqll. 
The number of diagrams contributing to is 36. The contribution to comes from 
252 diagrams (36(23 _ 1 ) = 252). 

Subset VI(b). Diagrams obtained by inserting a 114 in the internal photon lines of sixth-order 
diagram Mqh. The number of diagrams contributing to is 54. The contribution to 
comes from 162 diagrams (54(2^ — 1) = 162). 
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vi(j) VI(k) 


FIG. 6; Set VI 

Subset VI(c). Diagrams obtained by inserting a 112 in the photon lines connecting the closed 
lepton loop A 4 and the open lepton line in the eighth-order diagrams The number of 

diagrams contributing to is 144. The contribution to comes from 432 diagrams 
( 144(22 _ 1 ) = 432 ). 

Subset VI(d). This subset consists of diagrams in which the open lepton line of Mqh is 
modihed by fourth-order radiative corrections. The number of diagrams contributing to 
is 492. The contribution to comes also from 492 diagrams. 

Subset VI(e). This subset consists of diagrams in which the open lepton line of Mqh is 
modihed by second-order radiative correction whose photon line has 112 insertion. The 
number of diagrams contributing to is 48. The contribution to comes from 144 
diagrams (48(2^ — 1) = 144). 
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Subset VI(f). This subset is derived from diagrams of Msiib in which photon lines connecting 
Ag to the open lepton line receive a 112 insertion. The number of diagrams contributing to 
is 180. The contribution to comes from 540 diagrams (180(2^ — 1) = 540). 

Subset VI(g). This subset consists of diagrams in which both the closed loop A 4 and the 
open lepton line are modihed by second-order radiative corrections. The number of diagrams 
contributing to is 480. The contribution to comes also from 480 diagrams. 

Subset VI(h). This subset consists of diagrams in which A 4 is modihed by fourth-order 
radiative corrections. The nnmber of diagrams contributing to is 630. The contribntion 
to A 2 ^'^^ comes also from 630 diagrams. 

Subset VI(i). This subset consists of diagrams in which A 4 is modihed by second-order radia¬ 
tive correction whose photon has a v.p insertion 112 . The nnmber of diagrams contribnting 
to A^^^^ is 60. The contribntion to A^^^^ comes from 144 diagrams (60(2^ — 1) = 180). 

Subset VI(j). This snbset is obtained from Mqll (which contains a A 4 ) by hanging a second 
A 4 on two of three photon lines attached to the open lepton line. The nnmber of diagrams 
contribnting to A^^^^ is 54. The contribntion to comes from 162 diagrams (54(2^ — 1) = 
162). 

Subset VI(k). This subset consists of vertex diagrams containing Ag hve of whose photon 
lines end np on the open lepton line. The nnmber of diagrams contribnting to is 120. 
The contribution to comes also from 120 diagrams. 

The total nnmber of diagrams of Set VI contribnting to is 2298. The nnmber of 
diagrams of Set VI contributing to A^^^^ is 3594. 

The total nnmber of Feynman diagrams contribnting to A[^°^ is the sum of contributions 
from all diagrams described above, which is 12672. The total nnmber for A^^^^ is 9080. 

The nnmber of Feynman diagrams contribnting to Ag^*^^ can be readily derived from the 
above result. 

III. LEADING DIAGRAMS CONTRIBUTING TO 

Fortunately, it is not difficult to identify diagrams which may give large contribntion to 
They are diagrams containing ln(m^/me) terms which tend to have large nnmerical 
valnes because is much larger than rUe- 
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One source of ln(m^/me) is the vacuum-polarization contribution to the photon propa¬ 
gator which yields the logarithmic factor as a consequence of charge renormalization. The 
renormalized photon propagator has the form 




Da'ia) = £») + ■•■, 


(18) 


where, to order a, 


= 1 + 


a 


TT 


lln((/Vm^) - 5 + 


(19) 


When Dfi is inserted m. g — 2 diagrams, the momentum scale is set by the muon mass. Thus 
the a/vr term will give a factor of the order of 


2 5 

Kri = g ln(r((mp/m,)) - 


( 20 ) 


where rj is an fuzzy factor of order 1. — 3 for p = 1. 

As a matter of fact, more important sources of ln(m^/me) are diagrams built up from the 
large sixth-order diagram Mqli containing an external light-by-light scattering subdiagram 
of closed electron loop A 4 . (By external we mean that one of four photon leg is the external 
magnetic held. If all photons are virtual photons, we call it internal.) The extraordinary 
size of Mqii was initially discovered by numerical integration. The primary cause of large 
size is that it has a logarithmic mass-singularity for irie 0. But this is not the whole 
story. It was pointed out by Yelikhoskii Q that, in the large m^/me limit, the muon 
line may be regarded as a static source of Coulomb photons as well as hyperhne spin-spin 
interaction. Of three photons exchanged between the muon line and the electron loop, one 
photon is responsible for the spin-spin interaction while the other two are essentially static 
Coulomb potential. Integration over these Coulomb photon momenta gives a factor m each, 
contributing a factor 7 r^(~ 10 ) to the leading term 

aAY”'” = ^ Hmjm,). (21) 

For the physical value of this is about 35. Mqh as a whole is reduced by the 

negative mass-independent term to about 21, which is still very large. The value of Mqll 
containing an electron loop A 4 is known exactly and its numerical value is 


Mqll = 20.947 924 34 ( 21 ), 


( 22 ) 
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where the uncertainty is due to that of the muon mass only. 

The difference between Mqll and j^iay be interpreted as an indicator of the 

degree to which the picture of static Coulomb potential is valid. One way to incorporate 
this is to introduce a fudge factor ^ such that 


MeLL = 


(23) 


In this interpretation we have ^ ~ 0.77. 

These two sources of ln(m^/me), light-light-scattering loop and vacuum-polarization loop, 
can work together and give even larger numerical factors. For instance, the leading term of 
the integral A 2 \yI{a)\, which has insertion of two electron loops 112 and should have been 
written as VI(a)[e,e,e] following the notation of Table HVl will be of order 


A^lVKa)] ~ SA-pfsti, 


(24) 


where the factor 6 is the number of ways two electron loops 112 can be inserted in Mqll. 
This leads to A 2 [VI{a)] ~ 1130 for = 1 ~ 3). As is seen later, the actual value is 

about 543. Thus ~ 2. 

Similarly, the subset VI(b), which has insertion of one 114 in Mqll (VI(b)[e,e] in the 
notation of Table IHD, will give contribution of the order 

A2[VIib)] ~ 3 X ^ X K^Mqll ^ 142, (25) 

for 7] = 1 with the help of the identity 

n4('=^) = ^|n4P) + (^)h^(c(3) + T). (26) 

The numerical evaluation give A 2 [VI(b)]~ 169. Thus 77^ ~ 3.5 in this case. 

If we apply blindly the argument based on (PU|). AVG ^vou-lci obtfliin^ for 3^ 


A2[C/(c)]~ 3K,Msllc^27, 

A 2 [VI{e)] ~ K^Msllc^9, 

A2[VI{f)] ~ 3K^MsLLb ^ -4.5, 

A2[1"/(7)]~ K^MsLLb^-lA. (27) 


Unfortunately, these estimates are completely misleading because Mg^Lc is the sum of large 
terms which tend to cancel each other. Similarly for 
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Subsets VI(d), VI(g), and VI(j) have no 112 insertion and thus will have no particular 
enhancement although individual members of these subsets (being gauge-dependent) might 
have rather large values. 

Among the diagrams of Sets I - V, the term with the highest power of logarithm is 

^2[/(a)] ~ (28) 

although its value (~ 80) is not dominating because of its small numerical factor. Actually 
terms of lower logarithmic powers reduce this to an even smaller value (~ 20). 

In the study of eighth-order terms, the value of was found to be smaller than 3 in 
most cases and varies in the range 2 < < 2.5, except for VI(b). Thus the estimates given 

above for r] = 1 (or ~ 3) are likely to be overestimates by 20 ~ 50%. 

The contribution of the subset VI(k) had been estimated in a different manner. It is 
that the leading term in the large m^/rue limit is of the form 


A2\yi{k)\ = 7r'‘(0.438.. ln(m^/me) + ...). 


Based on this observation it was estimated that 



A2[VI{k)] ~ 185 ±85. 


(29) 


(30) 


This term is large mainly because of the presence of the vrd factor. Its origin can be readily 
understood by a generalization of the argument leading to Eq. (21) to the case in which 
2n±l photons are exchanged between the light-by-light type loop A 2 n+ 2 , n=l,2,..., and the 
muon line [44]. In the large m^/rrie limit this mechanism generates the structure 


^(2n+l) ^ -_ 


Numerical values of some c„ are known 



Cl 


C2 


2 

3’ 

0.438 .... 


(31) 


(32) 


IV. ANALYTICALLY KNOWN CONTRIBUTION TO a).'"’ 

At present only a small number of integrals in the subsets 1(a), 1(b), 1(c), 11(a), and 11(b) 
are known analytically. Their expansion in the ratio me/m^ are given in From Table 
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2 of 


4lj | we obtain 


a^[I{a)] = 22.566 973 (3), 
a^[I{b)] = 30.667 091 (3), 
a^[/(c)] = 5.141 395 (1), 

a^[II{a)] = -36.174 859 (2), 
a^[//(5)] = -23.462 173 (1), 


(33) 


where the uncertainties come from the measurement uncertainty of me/m^ only. These 
results show strong cancellation among diagrams of Sets I and II, which is analogous with 
the cancellation among Group I and Group II diagrams contributing to the eighth-order 
term of muon g-2. 


V. NUMERICAL EVALUATION OF DIAGRAMS CONTRIBUTING TO 

As was discussed in Section UTl diagrams containing light-by-light-scattering subdiagrams, 
namely the Set VI, is the source of dominant contribution to . However, let us describe 
the results of numerical evaluation starting from the Set I. 


A. Set I 


Integrals for the diagrams of Set I, except for subsets I(i) and I(j), can be readily ob¬ 
tained from eighth-order diagrams belonging to Group I by insertion of v-p loops. Diagrams 
numerically evaluated thus far are shown in Tables mun and 11111 Summing up the hrst four 
lines of Table H] one obtains 


A-Am = + Afgsr' + vgsr-”’ + MAjSi 


(e,e,e,m) 


r(e,e,m,m) 


r(e,m,m,m) 


= 22.567 05 (25). 

Similarly, we hnd from the rest of Table Q] 


Him 

^2 [/(C)] 


p4,p2:2 “T ^'^2,p4,p2:2 “C ^'^2,p4,p2:2 “r ^'^2,p4,p2:2 “r ^'^2,p4,p2: 

)67 54 (33), 


+ M 


-(e,m) 




5.141 38 (15), 


(34) 


(35) 
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The values of A 2 [I{a)\, A 2 [/(&)], and A 2 [I{c)] are in good agreement with semi-analytic 


values 


^1 quoted in 


A 2 [I{dy\ is evaluated from the entries in Table ITTl together with Table III and Table VI of 
1^ . which list residual renormalization terms. We obtain 


,Je(e),e) 

2,p4{p2)p2 


,^(m(e),e) 

2,pA{p2)p2 


^'^2,p4(p2)p2 


j.Am(m),e) 

2,pA{p2)p2 


o , a j.Ae{e),e) _ A \ R(e,e) j.Am,e,e) 

^^^''2,p4A(p2)p2 ^^^^2,p4B{p2)p2 ^^^2,p2 ^^-'2,p2:2 

7.45173 (101), 

9 A/r(™b)-e) I A AAm(e),e) _ A /\ R{m,e) 

^^'^2,p4A(p2)p2 ^^'^2,p4S(p2)p2 ^‘^-°2,p2 ^'"-'2,p2:2 

1.02576 (12), 

9 hAAm),e) A ,^(e(m),e) . , R(e,m) AAm,e,e) 

^^2,pAA{p2)p2 + 4T?2,p4S(p2)p2 “ ^^^2,p2 ^^2,p2-.2 

0.13084 (2), 

9 l\yr(mim),e) . A j.Am{m),e) _ A \ nim,m) AAm,m,e) 

^^^''2,p4A(p2)p2 ^^'^^2,piB{p2)p2 ^^^2,p2 ^''^2,p2:2 

7.17353 (81) X 10“^ 


(36) 


j.Ae(e),m) 

^'^2,p4(p2)p2 


AAm(e),m) 

^^■'2,p4(p2)p2 


AAe{m),m) 

^'^^2,p4(p2)p2 


9 i\Ae{e),rn) . j.Ae(e),m) a a p(e,e) j.Am,m,e) 

^^2,pAA{p2)p2 + 4tW2_p4B(p2)p2 “ ^^^2,p2 ^2,p2:2 


0.15845 (3), 

4.82662 (1) X 10-2, 

K 

5.29020 (63) x 10-^ 


9 _L A A/f(™h)-m) A A p{m,e) j.Am,m,m) 

^^2,pAA{p2)p2 + 4^2,p4S(p2)p2 “ 4AT?2,p2 ^2,p2:2 


9 A/fhM-™) I /I nAe{m),m) _ a /\ n(e,rn) AAm,e,m) 
^^'^2,p4A(p2)p2 + ^^'^2,p4S(p2)p2 ^^^2,p2 ^'"-'2,p2:2 


(37) 


From these contributions we obtain x42[/((i)]: 


4.[/(<i)i = a 4S7V+AfiSSi'h + + "bw 


(m(e),e) 


r(e(m),e) 




+ 


,^(e(e),m) ,^(m(e),m) ,Ae{m),m) 

^‘^^2,p4(p2)p2 ^‘^2,p4(p2)p2 ^'^2,p4(p2)p2 


= 8.89207 (102). 


(38) 


Several terms of x42[/((i)] was also evaluated by an alternative method using an exact spectral 
function of 114 ( 112 ) in which two loop masses are equal. They are listed at the bottom of 
Table HH They are in good agreement with the values given in (jSED and (jSZD. 

x42[/(e)] is evaluated using the Fade approximant for Ifg, which is known to be a very 
good approximation from an earlier work. Our result is 
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A2[I(e)] = h4% + 

= -1.219 20 (71). 


(39) 


^2[7(/)] is evaluated using the entries of Table HTTl together with Table IVTl We obtain 


(m,e(e,e)) 

®2,p4(p2:2) 


(m,e(e,m)) 

“2,p4(p2:2) 


(m,e(m,m)) 

®2,p4(p2:2) 


(m,m(e,e)) 

®2,p4(p2:2) 


®2,p4(p2:2) 


^^2,p4A{p2-.2) + ^^2,p4B(p2:2) “ ^^^2,p2-.2 ^2,p2 


^(m,e(e,e)) 


)(e(e,e)) j,Am,e) 


2.88598 (9), 

^^'^2,p4A(p2:2) 

0.16111 (3), 

^^'^2,p4A(p2:2) 

0.01063 (1), 

^^'^2,p4A(p2:2) 

0.53660 (9), 

^^'^2,p4A(p2:2) 

0.09079 (2), 


A A/r(™’®h’”^b 1 A _ 9 A A/fl*"’®) 

^^'^2,p4A(p2:2) + ^^'^2,p4B(p2:2) ^^^2,p2:2 ^^^2,p2 


A A/f(™’®(™-’™)) I A A/fl™’®!™-’™)) _ 9A D(e(«^.”^)) A/f(™>®) 
^^'^2,p4A(p2:2) ^^'^2,p4S(p2:2) ^^^2,p2:2 ^^^2,p2 


^^^2,p4A(p2:2) + ^^^2,p4B(p2:2) “ ^^^2,p2:2 ^^h,p2 


+ ^'^h,piB(p2-S) - ‘=“-t' 2 ,p 2,2 “ 2 ,p 2 


(40) 


Fi’ora these contributions we obtain Ac)/)/)]: 


42[/(/)] 


(m,e(e,e)) . (m,e(e,m)) . (m,e(m,m)) . (m.mfe.e)) . (m,m(e,m)) 

«gp4(p2:i + «2,p4(p2:2) + «gp4(p2:2) + «gp4(p2:2) + «2,p4(p2:2) 

3.68510 (13). 


(41) 


B. Set II 


Integrals of Set II, except for subset 11(e), can be readily obtained from eighth-order 
diagrams belonging to Group II by insertion of v-p loops. We have evaluated only subsets 
11(a), 11(b), and 11(f) thus far. From Table HVl and Table EH of this paper and Tables III 
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and VI of 


we obtain 


m: 


(e,e,e) _ r) ^Ae,e,e) 
4,p2:3 ~ ^^^''4a.o2:3 


M, 


(e,e)(e) 

4,p2:2,p2 


M, 


M, 


f ( 6 , 6 , 6 ) jyAe,e,e)a , ^Ae,e,e)l 3 a „ TyAe,e,e) 

4a,p2;3 + ^^46,p2:3 + ^^4b,p2:3 “ ^^2tW2,p2:3 

■28.43132 


(344), 


(e,e,m) 


4,p2:3 


(e, 6 )(m) 

4,p2:2,p2 


9/i//-(e,e)(e) i /i^(e>e)(e) , Ti^{e){e,e) _ A _ A A/f(®) 

^AV24„ p2:2,p2 “r ^‘^4fe,p2:2,p2 “r ^‘^4fe,p2,p2:2 ^^2,p2^^^2,p2:2 ^^2,p2:2^^^2,p2 

-27.42432 (127), 

9/l//-(e,6,m)() ,^(6,6,m)() ,^()(6,6,m) Ap ^Ae,e,m) a p(6,6,m) 

^^^ha,p2-.3 + J-^hb,p2:3 + ^^4b,p2:3 “ '^-02^142,p2:3 “ ^^2,p2:3 

9/l/f(e.e)(m) ,Ae,e)im) ,^(m)(6,e) _ \ n{i^) aA^^A _ A T\A^'> 

^^''^Aa,p2-.2,p2 ' ^^-'46,p2:2,p2 “r ^^-'46,p2,p2:2 ^■‘^2,p2^^-‘^2,p2:2 ^■‘^2,p2:2^^-‘-2,p2 


-6.79245 (56), 

'■^^^4a,p2-.2,p2 ^ 

-1.95703 (37), 


(42) 


A^(e,m)(6) 

^'^4,p2:2,p2 


M, 


.,^( 6 ,m)( 6 ) . c)j.Ae,m){e) , 9 ;.^( 6 )(e,m) 

^^^■'4a,p2:2,p2 “r ^^^■'4fe,p2:2,p2 ^^^-'4b,p2,p2:2 ' 


(e,m,m) 


4,p2:3 


j.Am,m){e) 

^^^4,p2:2,p2 


]\Ae,m){m) 

^^■'4,p2:2,p2 






= -4.14893(48), 


= 2M, 




4a,p2:3 


+ M, 


(e,m,m)() 


46,p2:3 


+ M, 


()(e,m,m) 


4fe,p2:3 


- A5oM 




2,p2:3 


- A5, 


(e,m,m) 

2,p2:3 


Mo 


= -0.95284 (14), 

r(m,m)(e) 


r) A/f('rn,'rrL)((i) , ,j(e)(m,m) _ A n(e) _ A j\/fy 

^^'^^Aa,p2-.2,p2 ' ^'^46,p2:2,p2 ^'^46,p2,p2:2 ^■‘^2,p2^'^^2,p2-.2 ^-'^2,p2:2 2*'22,p2 




?(e) l\Am,m) 


,(m,m) ,^(6) 


= -0.28902 (29), 

— 4 -L 9 _L 9 

— ^^^^Aa,p2-.2,p2 W Ab,p2-.2,p2 Ab,p2,p2:2 

= -0.47576(23). 


(43) 


Adding np these contribntions one obtains 

A2[//(a)] = -70.4717 (38). 


(44) 
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From Table El and Table EH of this paper and Tables III and VI of 



we obtain 


m: 


tplp2 — 2AiM4®’p4p2 + AiM4^p4p2 + + (p4 p2) 


M, 


(e,m) 


M, 


-(m,e) ^ . nj(rn,e)() . nj(rn,e)() . j^y)(m,e) / 4 

4,p4p2 ■^^l-^^-^4a,p4p2 ' Ab,pAp2 ' 4h,pAp2 ' \P^ V^j 


r(e.e)() 


(e,e)() 


r()(e,e) 




4,p4p2 


-19.04191 (203), 

+ A,Mi:-)? + + (p4 


AB^M, 


(e,m) 


^2,p4p2 

-1.28845 (25), 

r(m,e)() 




(m,e)() 


A)(m,e) 


p2) 


ABoM. 


(m,e) 


^2,p4p2 

= -3.13132 (25), 




(46) 


M. 


(e)(e) 


4,p4p2 


j.Ae){m) 

^^■'4,p4p2 


M, 


(m)(e) 


4,p4p2 


2A2A/iO]^p2 + ^2AfAp2 + (l>4 ^ p2) 

^bQM% - ab!,%a4% 

-9.42667 (145), 

2A2Afi:;^p'2 + A2A/^S2 + (P4 ” P2) 
AB < 2 a /^,!2 - ABgpA/a 

-0.60877 (18), 

2A2Afi”«'p>2 + A2A43;>2 + (p4 « p2) 

abIIa/S - 

-1.27436 (27). 


(46) 


Adding up these values we obtain 

A2[II{b)] = -34.7715 (26). (47) 

Table rvTTl lists numerical results of set 11(f) obtained in Version A. The results obtained 
by an alternate formulation (Version B) are listed in Table IVIIII Combining the values in 
these tables statistically one obtains 

A2[II{f)P^^ = -57.0633(109), 

As [//(/)] ("’"^^ = -4.7157(31), 

As[//(/)]= -15.6857( 37). (48) 
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From 


one obtains 


A2[IIif)] = -77.4648(120). 


(49) 


C. Set VI 

By far the largest contribution to comes from the subset VI(a), followed by the 

subset VI(b). Their integrals can be readily obtained by insertion of v-p functions 112 and 
n 4 into the sixth-order diagram Mqll- We have evaluated them precisely by VEGAS. The 
results are listed in Table HXl Summing them up, we obtain 


and 


A2[VI{a)\ = 629.1407 (118), 


A2[VI{h)] = 181.1285 (51). 


(50) 


(51) 


The difference between VI(a)[e,e,e] of Table HXl and old one in Eq.(2.49) of is due to a 
program error in the latter. The results in Table IIXI show that the leading term estimate 
(EH) is an overestimate by a factor two, while the estimate (ESI) is not too far off. 

Contributions of subsets VI(c), VI(e), VI(f) and VI(i) are also not difficult to evaluate, 
since they can be readily obtained by insertion of 112 into some eighth-order integrals. Their 
numerical values are listed in Tables |3 IXIl IXIIl and IXIIIl From this we obtain 


A2[EJ(c)](^’^) = 
A2[V/(c)](™’") = 
A2[VI{c)f'^^ = 
A2[VI{c)] = 


-17.0505 (1122), 

AA4:i% + Aa4TL% + AA4T,% + AA4’'4i% + AA4T,% - 

-14.2744 (105), 

+ A^M^lTIp + ^M^lTgp + + A^M^lTIp - 

-5.2514 (129), 

As [E/(c)] + A 2 [VJ(c)](”^’") + A 2 [VJ(c)]("’™) 

-36.5763 (1141), (52) 
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A2[VI{e)f^^^ 

A2[VI{e)f^^^ 

A 2 [VI{e)] 

A2[VlU)t^^^ - 

A2[VI{f)]^^^-^ : 
A2[VI{f)t’^'^ : 
A2[VIU)] ^ 

A2[VI{i)f’^^ = 
A2 [VI{i)]^^^^^ = 
A2[VI{i)]^^^^^ = 
A2[VI{^] = 


AA^rA + + AAf£f„ 


2AB'”J'm, 


6LL 


0.7524 (1338), 

AAC^4 + AA4’;i%^ + AAC41, + AACfi, + AAtdl 

2ABgi''A4”;”' 


-'6LL 

-3.9789 (78), 




- 2A4X^^6LL 
= -1.0950 (35), 

= A2[VI{e)f^^^ + A2[l^/(e)](™’") + A2[V 
= -4.3215 (1341), 

- 3AB2M^,l% 

-45.0425 (1463), 

- 3A52m(^;) 

8.7673 (228), 

\ I \]\ ^ L^,v, 

‘■SLLCp ^ ^^''■^SLLDp 


(53) 


AMifTL + AMifTL + AM4TL + AMifTL - 3ABM‘£Ti 


‘■SLLAp ^ ‘-^^'^^SLLBp 

-1.8847 (155), 

A2[VI{f)]^^’^^ + A2 [!"/(/)] + A2[VI{f)]^^’'^^ 

-38.1598 (1488), (54) 

- ‘iAB^l^M.Li 

28.3367 (1142), 

- 3ABi-fM^-r^ 

,0977 (105),, 

+ AMsS, + Am£S^ + Am£S, - SABil^^M,,, 

1.0983 (31), 

.2 [l//(z)] + A2[VI{t)]^^'^^ + A2[VI{t)]^^'"^^ 

27.3373 (1147). (55) 
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FIG. 7: Subset VI[j] 

The last columns of Table Hand Table IXTTl show that the factor r is roughly equal to SKrj 
with Kri ~ 2 for sets VI(c) and VI(f). These factors vary from diagram to diagram within 
each set, making naive estimates for the sums of diagrams given in (EZD entirely different 
from the corresponding terms of (jiH) - (I33j) obtained by explicit numerical integration. The 
’’enhancement” factor is r'U 1 for sets VI(e) and VI(i). Table IXTl and Table IXlTTl show 
that Krj ~ 1 in these cases, indicating that there is no signihcant enhancement. This is 
probably because 112 is buried deep in other subdiagrams in these cases. 

Among the remaining subsets the subset VI(k) is likely to be the largest and most chal¬ 
lenging because it has a subdiagram Ag which consists of 120 proper lepton loops to which 
six photon lines are attached. See Fig. 8. This is a diagram which appears for the hrst time 
in the tenth-order so that it cannot be derived from or related to lower-order diagrams. 

Another point of interest is that its leading logarithmic term is known PI|) and was used 
to estimate its contribution to However, to determine the actual contribution of VI(k) 
we must find the nonleading term, too. To answer this question it is best to evaluate VI(k) 
explicitly. 

It turned out that this is not difficult. Our first step is to reduce the number of inde¬ 
pendent integrals to 12 using the Ward-Takahashi identity, and reduce it further to 9 using 
the time-reversal symmetry pll|. Each integral generated by FORM has more than 90,000 
terms occupying about 30,000 lines of FORTRAN code. This is certainly huge, but not 
unmanageable since it is only 30 times larger than typical eighth-order integrals. We also 
note that this subset is particularly simple in the sense that it is entirely free from UV- and 
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FIG. 8: Subset VI[k] 

IR-divergences. Numerical integration over 13-dimensional Feynman parameter space can 
be handled without difficulty by VEGAS. The result of numerical integration is listed in 
Table 1X3 from which we obtain 


yl2[V/(fc)] = 97.123 (62). (56) 

Clearly the previous estimate was an overestimate by about 100. 

Another possibly large term is VI(j) [162 vertex diagrams] which we decided to evaluate 
explicitly. See Fig. 7. With the help of Ward-Takahashi identity and time-reversal invariance 
it can be represented by four independent integrals. FORM generated about 42,000 terms 
for each integral occupying about 18,000 lines of FORTRAN code. The result of numerical 
integration is listed in Table IM from which we obtain 


A2[VI{3)] = -25.505 (20). 


(57) 


The subsets yet to be evaluated are VI(d), VI(g), and VI(h). We foresee no technical 
problem in dealing with these subsets. 

Some results described in this section can be compared with the results obtained by the 
renormalization group method. See y for details. 
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VI. CONTRIBUTION TO ELECTRON G-2 


All Tables also contain values of mass-independent contributions from 958 vertex dia¬ 
grams belonging to 17 gauge-invariant subsets. These are actually contributions to 
namely the electron g-2. Their values, including residual renormalization terms whenever 


they are required, are listed below: 

Ai[J(a)] = 
A,[I{b)] = 

MHc)] = 
A,[I{d)] = 

Mm = 
MHf)] = 

Til [//(a)] = 
A,[II{b)] = 

MI I if)] = 

amiM 

A,[VI{b)] 

Ai[VI{c)] 

A,[VI{e)] 

Ai[VI{f)] 

Ai[Vm] 

Ai[VI{j)] 

Ai[VI{k)] 

VII. DISCUSSION 


4.7094 (6) X 10■^ 

7.0108 (7) X 10■^ 

2.3468 (2) X 10"^ 

4.4517 (5) X 10-^ 

1.0296 (4) X 10“^ 

8.4459 (14) X 10 -\ (58) 

4.130 (90) X 10-^ 

-5.422 (4) X 10“^ 

-2.434 (2). (59) 

= 1.0417 (4), 

= 1.3473 (3), 

= -2.5922 (34), 

= -0.4312 (6), 

= 0.7703 (24), 

= -0.0438 (11), 

= -0.2288 (17), 

= 0.6802 (38). (60) 


Let us first add up all terms contributing to A'^^\m^/me) of muon g-2 evaluated by 
numerical integration. From (jHll), (IH^ . dSHl), (EH), (EH), (EH, EH, (EH, (EH, (EH, (EH, 


(ESI), (ISH), (EH), (EH, EH, which represent the contributions of 2958 Feynman diagrams 
belonging to 17 gauge-invariant sets, we obtain 

me) [partial sum] = 662.50 (27). (61) 

The uncertainties from these diagrams, which include all dominant sources of uncertainties 
considered previously, have been reduced to an insignihcant level. Note, however, that some 
terms which were not included in previous estimates turned out to be not negligible. 

Of course, the real value of A^ 2 ^\mn/me) is not known until remaining diagrams are 
evaluated. However, they have no known mechanism for giving rise to large values and 
likely to remain modest in size and uncertainty. We therefore expect that the hnal value 
will stay within the range 

A^2^\mn/me)[estimate] = 663 (20). (62) 

This will reduce (fTH by 1.81 x 10“^^. 

Our next step is to evaluate 6122 vertex diagrams from the remaining 14 gauge-invariant 
subsets. Many of these diagrams can be integrated by means of available information on H 
and A. They include: 

1. Subsets 1(f), 1(g), 1(h) of Set I. 

2. Subsets 11(c), 11(d) of Set II. 

3. Subsets Ill(a) and Ill(b) of Set III. 

4. All diagrams of Set IV, which can be evaluated by simple modihcation of codes of 
Group V of eighth-order diagrams. 

These subsets are on our next time schedule. 

The remainder of diagrams are more difficult to evaluate for the following reasons: 

1. Subsets I(i) and I(j) require the knowledge of the eighth-order v-p spectral function 
Hg which have not yet been constructed. 

2. Diagrams of set 11(e) require construction of radiatively-corrected light-by-light scat¬ 
tering subdiagrams, which are not yet available. 
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3. Subset III[c] contains A 4 internally. 

4. Subsets VI(d), VI(g), and VI(li) have no lower-order structure upon which they can 
be built. 

We foresee no intractable barrier for their evaluation. 

As far as the electron g -2 is concerned what is most important is the Set V which have 
no lower-order fermion loop structure. The integrands of this set are gigantic and require an 
enormous number of UV and IR subtraction terms. Our experience with the subsets VI(k) 
and VI(j) indicates, however, that their sizes are still manageable with available computers. 
What is really crucial for their evaluation, however, is that all steps of construction of 
integrand must be fully automated. Thus far we have succeeded in obtaining a code that 
handles renormalization of ultraviolet divergence automatically For the moment the 
infrared divergence is treated by giving a small cutoff mass A to the photon. A code for 
cutoff-independent treatment of IR divergence is being developed. 
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TABLE I: Numerical evaluation of diagrams of subsets (a), (b), (c), and (e) of Set I. The notation 


follows that of 


42|| with some modification and adaptation, nj? is the number of Feynman diagrams 


represented by the integral. 


Integral 

Up 

Value (error) 

Sampling per 

No. of 



including rej? 

iteration 

iterations 

^^^2,p2:4 

1 

20.14293 (23) 

1 X 10® 

180 

f.Ae,e,e,m) 

^'^2,p2:4 

4 

2.20327 ( 9) 

1 X 10^ 

80 

^^^2,p2-A 

6 

0.20697 ( 2) 

1 X 10^ 

20 

^^2,p2:4 

4 

0.01388 ( 1) 

1 X 10^ 

20 

^^2,p2:4 

1 

4.7094 (6)xl0-^ 

1 X 10^ 

20 

^'^2,p4,p2:2 

9 

27.69038 (30) 

4 X 10® 

190 

, ,(e,m,e) 
^'^2,p4,p2:2 

18 

1.16628 ( 9) 

1 X 10^ 

100 

i\Ae,m,m) 
^'^2,p4,p2:2 
, .(m,e,e) 
^'^2,p4,p2:2 

9 

0.03182 ( 3) 

1 X 10® 

20 

9 

1.61436 ( 6) 

1 X 10^ 

120 

^'^2,p4,p2;2 

18 

0.16470 ( 5) 

1 X 10® 

20 

■'^2,p4,p2:2 

9 

7.0108 (7) X 10-3 

1 X 10^ 

20 

^'^2,p4:2 

9 

4.74212 (14) 

1 X 10® 

220 

^'^2,p4:2 

18 

0.39926 ( 3) 

1 X 10^ 

120 

,Am,m) 

^'^2,p4:2 

9 

2.3468 (2) X 10-2 

1 X 10^ 

20 

^'^2,p6p2 

30 

- 1.20841 (70) 

1 X 10® 

100 

^'^2,p6p2 

^'^2,p6p2 

30 

-0.02110 ( 4) 

1 X 10^ 

120 

30 

0.01031 ( 1) 

1 X 10® 

20 

,Am,m) 

^'^2,p6p2 

30 

1.0296 (4) X 10-2 

1 X 10^ 

20 


33 








TABLE II: Numerical evaluation of diagrams of subsets (d) of Set I. The notation follows that of 
with some modification and adaptation, np is the number of Feynman diagrams represented 
by the integral. Subscripts pAA, pAB refer to part of fourth-order u — p 114 containing vertex 
correction and self-energy insertion, respectively. The last four lines are values obtained using the 
exact sixth-order spectral function 114(2) • 


Integral 

Tip 

Value (error) 

Sampling per 

No. of 



including np 

iteration 

iterations 

Ji/r(e(e).e) 

2,pAA{p2)p2 

2 

2.63064 (72) 

1 X 

10^ 

80 

. .(m(e),e) 
^'^2,p4U(p2)p2 

2 

0.76997 ( 7) 

1 X 

10^ 

40 

. ,(e(m),e) 
^’^2,p44l(p2)p2 

2 

0.12703 ( 2) 

1 X 

10^ 

40 

^'^^2,p4:A{p2)p2 
. ,(e(e),m) 
^'^^2,p4A{p2)p2 

2 

7.3352 ( 7) X 10-2 

1 X 

10^ 

40 

2 

1.3211 (18) X 10-2 

1 X 

10^ 

79 

2,pAA{p2)p2 

2 

3.5428 ( 4) X 10-2 

1 X 

10^ 

40 

2,p4:A(p2)p2 

2 

5.0276 ( 6) X 10-3 

1 X 

10^ 

40 

^'^2,p4A(p2)p2 

2 

3.8330 ( 4) X 10-3 

1 X 

10^ 

40 

^r(e(e),e) 

^'^^2,p4B{p2)p2 

4 

5.51053 (70) 

1 X 

10^ 

80 

j.Jm{e),e) 

^’^2,p4S(p2)p2 

4 

0.63490 ( 9) 

1 X 

10^ 

40 

,Ae{m),e) 

^'^2,p4B(p2)p2 

4 

3.9114 (47) X 10-3 

1 X 

10^ 

29 

2,p4B(p2)p2 
. ,(e(e),m) 
^’^2,p4S(p2)p2 

4 

1.1129 ( 4) X 10-2 

1 X 

10^ 

29 

4 

0.15798 ( 2) 

1 X 

10^ 

80 

2,p4S(p2)p2 

4 

3.2133 ( 5) X 10-2 

1 X 

10^ 

40 

^'^2,p4B(p2)p2 

4 

2.6449 (19) X 10-^ 

1 X 

10^ 

40 

^’^2,p4S(p2)p2 

4 

6.1873 (17) X 10-^ 

1 X 

10^ 

40 

,Je(e),e) 

2,p4(p2)p2 

6 

7.45270 (88) 

1 X 

10® 

140 

,Je{e),m) 

2,p4(p2)p2 

6 

0.15853 (10) 

1 X 

10^ 

120 

2,p4(p2)p2 

6 

0.07173 ( 3) 

1 X 

10® 

20 

2,p4(p2)p2 

6 

3.8028 (5) X 10-3 

1 X 

10^ 

20 
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TABLE III: Numerical evaluation of diagrams of subsets (f) of Set I. The notation follows that of 
with some modification and adaptation, np is the number of Feynman diagrams represented 
by the integral. Subscripts pAA, pAB refer to part of fourth-order u — p 114 containing vertex 
correction and self-energy insertion, respectively. 


Integral 

np 

Value (error) 

Sampling per 

No. of 



including np 

iteration 

iterations 

.^(e(e,e)) 

^^^■'2,p4A(p2:2) 

I 

1.99747 (10) 

4 X 10^ 

240 

^^'^2,p4A(p2:2) 

2 

0.15863 (3) 

1 X 10^ 

20 

^^'^2,p4A(p2:2) 

1 

1.0612 (2) X 10"2 

1 X lO'^ 

20 

^^'^2,p4A(p2:2) 

^^'^2,p4A(p2:2) 

1 

0.44361 (6) 

1 X 10^ 

20 

2 

9.4542 (13) X 10"2 

1 X 10^ 

20 

^^^h,p4A{p2-.2) 

1 

8.1763 (13) X 10-3 

1 X 10^ 

20 

^^^^2,p4B{p2:2) 

2 

0.94959 (11) 

1 X 10^ 

340 

^^'^2,p4S(p2:2) 

4 

2.5573 (56) x 10-3 

1 X 10^ 

20 

A AAA'm,m)) 
^^'^2,p4S(p2:2) 

^^'^2,p4B(p2:2) 

2 

1.6448 (275) x lO-® 

1 X 10'^ 

20 

2 

0.26023 (6) 

1 X 10^ 

20 

^^^^2,p4B(p2:2) 

4 

1.1060 (5) X 10-2 

1 X 10^ 

20 

^^'^2,p4S(p2:2) 

2 

2.6957 (32) x 10-^ 

1 X 10^ 

20 
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TABLE IV: Numerical evaluation of set 11(a) which consists of Il(ai) and 11 ( 02 ). The set Il(ai) 
has contributions denoted by Ai and set 11 ( 02 ) has contributions denoted by A 2 . The symbol ()() 
indicates photon lines in which v-p loops 11 are inserted, np is the number of Feynman diagrams 
represented by the integral, t.r. is time-reversed Feynman diagram. 


Integral 

np 

Value (error) 

Sampling per 

No. of 



including np 

iteration 

iterations 

+ t.i'. 

6 

15.15611 (283) 

1 X 10® 

160 

A.MiiS’" + t-r- 

18 

1.31796 ( 49) 

1 X 10® 

200 

+ t.r. 

18 

0.09109 ( 13) 

1 X 10® 

100 

AiACS"”’ + ‘'t- 

6 

0.00413 ( 9) 

1 X 10^ 

20 

AlAfiSrf’ +AiAA«;«' 

6 

-30.09130 (191) 

4 X 10 '^ 

200 

^lT'4fe,p2:3 +^l^'^46,p2:3 

18 

-6.38188 ( 26) 

1 X 10® 

100 

^lT'4fe,p2:3 +^l^'^46,p2:3 

18 

-0.84226 ( 4) 

1 X 10® 

120 

AiA<™r>‘’ +AiAa“23-'’”' 

6 

-0.05422 ( 4) 

1 X 10 '^ 

20 

A 2 Mi:;;gp 2 + t.r. 

6 

12.23457 (99) 

2 X 10® 

370 

^2^'o4a,p2:2,p2 + 

6 

0.01971 (34) 

1 X 10^ 

120 

^ 221040 ,p 2 : 2 ,p 2 + 

12 

0.27276(35) 

2 X 10 '^ 

380 

A A/r(™’™he) , r. „ 

^22'o4a,p2:2,p2 + 

6 

-0.01408(27) 

1 X 10 '^ 

20 

^22W4a^p2:2,p2 + 

12 

-0.08222(22) 

1 X 10^ 

20 

^ 221040 ,p 2 : 2 ,p 2 + 

6 

-0.00998( 3) 

1 X 10 '^ 

20 

^^^Ab,p2:2,p2 +^^^hb,p2,p2:2 

6 

-28.69942 ( 79) 

1 X 10® 

380 

^‘i^Ab,p2-.2,p2 +^‘i^Ab,p2,p2-.2 

6 

-1.72076 (14) 

1 X 10 '^ 

100 

A _i_A AAAie,m) 

^2^Ab,p2:2,p2 +^‘2^Ab,p2,p2-.2 

12 

-3.71561 (32) 

1 X 10 '^ 

120 

^2^Ab,p2:2,p2 +^^^hb,p2,p2:2 

6 

-0.23956 ( 8) 

1 X 10 '^ 

20 

^^^Ab,p2:2,p2 +^‘2^Ab,p2,p2-.2 

12 

-0.37976 ( 7) 

1 X 10 '^ 

20 

^‘i^Ab,p2:2,p2 +^‘i^Ab,p2,p2-.2 

6 

-0.03619 ( 1) 

1 X 10 '^ 

20 
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TABLE V: Numerical evaluation of diagrams of set n(b). Ai indicates contributions in which II’s 
act on the same photon line while A 2 refers to those in which IT’s act on different photon lines. 
The symbol ()() indicates photon lines in which v-p loops 11 are inserted, np is the number of 
Feynman diagrams represented by the integral, t.r. is time-reversed Feynman diagram. 


Integral 


np 

Value (error) 

Sampling per 

No. of 




including np 

iteration 

iterations 

^iMiZ £2 + iP^ ^ p2)+t.r. 


36 

12.65000 (135) 

4 X 10® 

220 

^ p2)+t.r. 


36 

0.30749 ( 22) 

4 X 10^ 

140 

^ p2)+t.r. 


36 

0.84809 ( 23) 

1 X 10® 

160 

^ p2)+t.r. 


36 

0.04424 ( 33) 

1 X 10^ 

20 

^iM^b,p}p2 +^^^2,ptp2 +(P^ ^ 

p2) 

36 

-19.70781 (143) 

1 X 10® 

220 

AiAC“ +A.MiK +(f-4 „ 

p2) 

36 

-1.23984 ( 10) 

4 X 10^ 

140 

AiA/i;y» +a.m“;“' +(„4 „ 

p2) 

36 

-3.06974 ( 9) 

1 X 10® 

200 

AiM^-'O +A.a4«5;4 +(j,4„j,2) 

36 

-0.30984 ( 15) 

1 X 10^ 

20 



18 

5.89471 (107) 

1 X 10® 

240 



18 

0.09286 ( 16) 

1 X 10^ 

140 



18 

0.21886 ( 23) 

1 X 10^ 

180 



18 

-0.01723 ( 5) 

1 X 10^ 

20 

+{pA^p2) 


18 

-9.83570 (97) 

4 X 10^ 

260 



18 

-0.56867 ( 8) 

1 X 10^ 

120 

^2Mi'^l^,+{p4^p2) 


18 

-1.36297 (13) 

1 X 10^ 

100 

+(p4^p2) 


18 

-0.10756 ( 2) 

1 X 10^ 

20 
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TABLE VI: Numerical evaluation of diagrams for auxiliary quantities needed to evaluate contri¬ 
bution of Set I and Set II. 


Integral 

np 

Value (error) 

Sampling per 

No. of 



including np 

iteration 

iterations 

A 

^'^2,2:3 

1 

16.15765 (75) 

1 X 10^ 

100 

A D(e,e,m) 

^-°2,2:3 

3 

2.71594 (12) 

1 X 10^ 

100 

^-°2,2:3 

3 

0.36136 ( 4) 

1 X 10^ 

40 

A riiTn^m^m) 

^-^2,2:3 

1 

0.02381 ( 1) 

1 X 10^ 

20 

A 

^-°2,p4p2 

6 

13.27621 (93) 

1 X 10^ 

100 

A 

^-°2,p4p2 

A r(™’®) 
^-°2,p4p2 

^-°2,p4p2 

6 

6 

6 

0.53279 ( 5) 

1.31884 ( 5) 

0.13066 ( 2) 

1 X 10^ 

1 X 10^ 

1 X 10^ 

80 

120 

20 

A R(e(e,e)) 

^-°2,p2:2 

1 

0.02791 ( 1) 

1 X 10® 

40 

A oim(e,e)) 

^-°2,p2:2 

1 

5.33035 (20) 

1 X 10^ 

100 

^-°2,p2;2 

2 

0.47208 ( 9) 

1 X 10® 

40 

A o(e(m,e)) 
^-°2,p2;2 

2 

3.6159 (8) X 10-5 

1 X 10® 

40 

A D(e(m,m)) 

^-°2,p2;2 

1 

8.3020 (36) X 10-^ 

1 X 10® 

20 
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TABLE VII: Numerical evaluation of diagrams of set 11(f) in Version A. The suffix p indicates 
insertion of 112 in the photon lines connecting the muon line and the I — I loop. = 

= AMsLLxp/AMsLLx, where x = J,K, L, in column 6 is for 
comparison with the enhancement factor for set 11(f) discussed in Sec. IIIII The logarithmic 
enhancement comes from the v — p loop only. This is consistent with r ~ 11 ~ 15 for (e,e), (me) 
and r ~ 0.7 ~ 2 for (e,m), (m,m). 


Integral 

np 

Value (error) 

Sampling per 

No. of 

r 



including np 

iteration 

iterations 



24 

70.6567 (254) 

4 X 10^ 

220 

11.05 


24 

35.0760 ( 78) 

2 X 10^ 

200 

13.74 


24 

6.0056 ( 26) 

2 X 10^ 

200 

0.94 

AMlSl 

24 

3.7717 ( 17) 

1 X 10^ 

140 

1.48 

^■^^^SLLKp 

24 

-87.8367 (291) 

8 X 10^ 

210 

11.29 

^^^^SLLKp 

24 

-26.1793 ( 83) 

2 X 10^ 

220 

13.97 

‘^^^^SLLKp 

24 

-5.4760 ( 29) 

2 X 10^ 

230 

0.70 


24 

-2.9311 ( 17) 

1 X 10^ 

140 

1.56 

^a4iL 

24 

-39.9514 (287) 

4 X 10^ 

310 

13.07 

AA/X 

24 

-24.5877 ( 81) 

2 X 10^ 

220 

14.73 

AA/£ri 

24 

-5.2468 ( 29) 

2 X 10^ 

200 

1.71 


24 

-3.2771 ( 17) 

1 X 10^ 

140 

1.96 


A2[2/](^’^) 

72 

-57.1314 (481) 

^2[2/]("A,e) 

72 

-15.6910 (140) 

A2[2/](^’™) 

72 

-4.7172 ( 49) 

Aa [2/] 

72 

-2.4365 ( 29) 
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TABLE VIII: Numerical evaluation of diagrams of set 11(f) in Version B. The symbol p indicates 
insertion of 112 in the photon lines connecting the muon line and the light —light loop. JKLp^^’y^ = 
AMg^’^jp + + AAfg^’^^p. The suffix 2 or 13 below indicates the muon line into which 

the magnetic field vertex is inserted. A2[2/](*’^i = JKLp^^'^^ + JKLp^^^^\ r = JKLp/JKL in 
column 6 is for comparison with the enhancement factor for set 11(f) discussed in Sec. II I II 


Integral 

Uf 

Value (error) 


Sampling per 

No. of 

r 



including np 



iteration 

iterations 


JKLp^^’^'> 

24 

-11.37774 (313) 

1 X 

10^ 

,1 X 10^ 

250, 20 

12.35 


24 

-1.87116 ( 30) 



1 X 10^ 

250 

15.86 


24 

-0.74946 ( 64) 



1 X 10^ 

250 

0.81 


24 

-0.28035 ( 39) 



1 X 10^ 

250 

2.38 

JKLp^^f 

48 

-45.68188 (1079) 

1 X 

10®; 

,1 X 10^ 

450, 100 

13.03 


48 

-13.81420 ( 386) 

1 X 

10^: 

,1 X 10^ 

250, 40 

15.82 


48 

-3.96527 ( 399) 



1 X 10^ 

250 

1.13 


48 

-2.15339 ( 182) 



1 X 10^ 

250 

2.47 


A2[2/](->-) 

72 

-57.0596 (113) 

A2[2/](”^’") 

72 

-15.6854 (39) 

A2[2/](^’”^) 

72 

-4.7147 ( 40) 

^2[2/]K™) 

72 

-2.4338 ( 19) 


40 







TABLE 
that of 


IX: Numerical evaluation of diagrams of Set VI(a) and Set VI(b). 
with some modification and adaptation. 


The notation follows 


Integral 

Tip 

Value (error) 

Sampling per 

No. of 



including np 

iteration 

iterations 

^'^^6LL,p2p2 

6 

542.91180 (910) 

1 X 10^ 

300 

j\Ae,e,m) 

^'^^6LL,p2p2 

72 

39.00349 (272) 

1 X 10® 

200 

]\Ae,m,m) 

^'^^6LL,p2p2 

36 

2.43029 ( 22) 

1 X 10® 

200 

. .(m,e,e) 
^^^6LL,p2p2 

36 

34.42389 (680) 

1 X 10® 

200 

j.Am,e,m) 

^^^6LL,p2p2 

72 

10.37125 (162) 

1 X 10® 

200 

,Am,m,m) 

^'^^6LL,p2p2 

36 

1.04171 ( 37) 

1 X 10® 

20 

^^^6LL,p4 

54 

168.72855 (478) 

1 X 10^ 

200 

^^^6LL,p4 

54 

7.58383 ( 50) 

1 X 10® 

210 

^^^6LL,p4 

54 

4.81614 (164) 

1 X 10® 

200 

^^^6LL,p4 

54 

1.34726 ( 30) 

1 X 10® 

140 
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TABLE X: Numerical evaluation of diagrams of Set VI(c). The notation follows that of m with 
some modification and adaptation. The suffix p indicates insertion of 112 in the photon lines 
connecting the muon line and l-l loop, r = AMsLLxp /where x = E,F, G, H, I, in column 
6 is for comparison with the crude enhancement factor for set VI(c) discussed in Sec. IIIII 


Integral 

Up 

Value (error) 

Sampling per 

No. of 

r 



including np 

iteration 

iterations 



18 

-82.9940 (141) 

2 X 10® 

200 

3.84 

<1 

18 

-1.8277 ( 18) 

1 X lO’^ 

40 



18 

-3.3463 ( 44) 

1 X 10^ 

140 



18 

-0.6168 ( 6) 

1 X 10^ 

140 


A ^(e.e) 
^■^^^SLLFp 

36 

-322.4493 ( 573) 

2 X 10® 

280 

4.25 

^■^^^SLLFp 

36 

-5.2571 ( 60) 

1 X 10^ 

80 


^^^^SLLFp 

36 

-9.3199 ( 64) 

4 X 10^ 

390 


^^^^SLLFp 

36 

-1.5392 ( 17) 

1 X 10^ 

280 



36 

-181.5345 (489) 

2 X 10® 

210 

5.17 

AMferL 

36 

-3.3841 ( 61) 

1 X 10^ 

80 


am‘S,'L 

36 

-6.5157 ( 84) 

4 X 10^ 

300 



36 

-0.9499 ( 18) 

1 X 10^ 

300 


^M^Tlhp 

18 

230.3344 (590) 

4 X 10® 

400 

4.27 

AMferL 

18 

2.2093 ( 74) 

1 X 10^ 

66 


AAfc'L 

18 

2.1519 ( 96) 

4 X lO’^ 

480 


AmSp 

18 

0.2839 ( 18) 

1 X lO’^ 

320 


A ^(e,e) 

36 

514.7317 (567) 

2 X 10® 

300 

4.57 

A^(e,m) 

^^^^SLLIp 

36 

7.0544 ( 59) 

1 X 10^ 

80 


A^(»n,e) 

36 

9.2489 ( 68) 

4 X 10^ 

320 


^■^''^SLLIp 

36 

1.1912 ( 13) 

1 X 10^ 

340 
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TABLE XI: Numerical evaluation of diagrams of Set VI(e). The notation follows that of m with 
some modification and adaptation. The suffix q refers to insertion of 112 in radiative corrections 
to the muon line, r = AMg^ixp /where x = E, F,G, H, I, in column 6 is for comparison 
with the crude enhancement factor for set VI(e) discussed in Sec. II I II 


Integral 

np 

Value (error) 

Sampling per 

No. of 

r 



including np 

iteration 

iterations 



6 

-35.1438 ( 514) 

1 X 10^ 

80 

1.62 


6 

0.0505 (16) 

1 X 10^ 

40 



6 

-0.9986 (17) 

1 X 10^ 

120 



6 

-0.1347 (2) 

1 X 10^ 

80 



12 

-100.5201 (458) 

4 X 10^ 

220 

1.33 


12 

0.3570 (19) 

1 X 10^ 

100 



12 

-2.4368 (41) 

1 X 10^ 

140 



12 

-0.2254 (3) 

1 X 10^ 

100 



12 

-38.1520 ( 399) 

1 X 10^ 

80 

1.09 


12 

-0.2662 ( 7) 

1 X 10^ 

40 



12 

-1.6144 (33) 

1 X 10^ 

180 



12 

-0.0957 (3) 

1 X 10'^ 

80 



6 

64.5209 ( 879) 

1 X 10^ 

80 

1.19 


6 

-0.4527 (12) 

1 X 10'^ 

40 



6 

0.3954 (39) 

1 X 10^ 

160 


aa4S', 

6 

-0.221 (91) X 10-3 

1 X 10^ 

80 



12 

189.0518 (620) 

4 X 10^ 

200 

1.68 

aa43 

12 

1.8726 ( 19) 

1 X 10^ 

80 


AAfe?, 

12 

2.0747 ( 37) 

1 X 10^ 

200 



12 

0.0719 ( 2) 

1 X 10^ 

80 
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TABLE XII: Numerical evaluation of diagrams of Set VI(f). The notation follows that of m 
with some modification and adaptation. The suffix p indicates insertion of 112 in the photon lines 
connecting the muon line and the I — I loop, r = AM^i^xp /where x = A,B,C,D, in 
column 6 is for comparison with the enhancement factor 3K^ for set VI(f) discussed in Sec. II I II 


Integral 

np 

Value (error) 

Sampling per 

No. of 

r 



including np 

iteration 

iterations 


A 

30 

307.3206 (848) 

4 X 10^ 

200 

5.90 

A 

30 

8.9175 ( 82) 

1 X 10^ 

120 


A 

^■^'^^SLLAp 

30 

2.9097 (118) 

1 X 10^ 

100 


^^^^SLLAp 

30 

0.7878 ( 10) 

1 X 10^ 

200 


<! 

60 

-482.5729 (603) 

1 X 10® 

200 

6.43 

AmIlSp 

60 

-16.0636 ( 78) 

1 X 10^ 

120 



60 

-13.7242 (121) 

2 X 10^ 

100 



60 

-1.7449 ( 12) 

1 X 10^ 

280 


<! 

60 

645.3472 (823) 

1 X 10® 

200 

6.00 

AMi“L 

60 

19.7833 ( 82) 

1 X 10^ 

180 


AAfc'L 

60 

30.4954 (125) 

2 X 10^ 

240 



60 

3.4824 ( 16) 

1 X 10^ 

300 


^^4Tld, 

30 

-252.4292 (616) 

4 X 10^ 

200 

6.67 

^^^^SLLDp 

30 

-8.4527 ( 66) 

1 X 10^ 

80 


^^'^^SLLDp 

30 

-1.1736 ( 86) 

1 X 10^ 

140 


^^^^SLLDp 

30 

-0.4074 ( 8) 

1 X 10^ 

200 
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TABLE XIII: Numerical evaluation of diagrams of Set VI(i). The notation follows that of m with 
some modification and adaptation. The suffix q refers to insertion of 112 in radiative corrections to 
the I — I loop, r = where x = A,B, C, D, in column 6 is for comparison with 

the enhancement factor for set VI(i) discussed in Sec. IIIII 


Integral 

Up 

Value (error) 

Sampling per 

No. of 

r 



including np 

iteration 

iterations 



10 

29.6915 (199) 

1 X 10^ 

80 

0.57 


10 

0.0914 ( 3) 

1 X 10^ 

40 



10 

0.4867 (30) 

1 X 10^ 

100 



10 

0.8068 (93) X 10-2 

1 X 10'^ 

140 



20 

-81.3369 (946) 

1 X lO'^ 

95* 

1.08 


20 

-1.2816 (26) 

1 X 10^ 

100 


AA&l 

20 

-4.6043 (69) 

1 X 10^ 

138 


Avis’, 

20 

-0.4926 (8) 

1 X 10^ 

200 


<1 

20 

71.0664 (409) 

1 X 10'^ 

80 

0.66 

aa4S, 

20 

0.7624 ( 7) 

1 X lO'^ 

40 


AMiS, 

20 

8.5729 (59) 

1 X lO'^ 

160 


AMiSl, 

20 

0.5664 (4) 

1 X 10^ 

100 


kC; 
oj' kC; 

X?® 

<1 

10 

-43.7735 (449) 

1 X 10^ 

80 

1.16 

AM<S, 

10 

-0.6699 (14) 

1 X 10^ 

60 


AMiS, 

10 

-0.2587 (42) 

1 X 10'^ 

100 


AMiS, 

10 

-0.0551 (5) 

1 X 10^ 

140 
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TABLE XIV: Numerical evaluation of diagrams of Set VI(j) contributing to the muon g-2. In 
the superscript (a, 6), a and b refer to the external and internal light-by-light-scattering loop, 
respectively. Last 4 lines with superscript (m, m) are mass-independent so that they are identical 
with the contributions to the electron g-2. 


Integral 

np 

Value (Error) 

Sampling per 

No. of 



including np 

iteration 

iterations 


24 

0.57928 (0.01386 ) 

1 X 10®,1 X 10^ 

100, 294 

XQjM 

12 

-16.91235 (0.00630 ) 

1 X 10®,1 X 10^ 

100, 100 


12 

-23.00801 (0.00777 ) 

1 X 10®,1 X 10^ 

100, 160 


6 

15.38181 (0.00388 ) 

1 X 10®,1 X 10^ 

100, 87 


24 

4.75211 (0.00609 ) 

1 X 10® 

105 

.(e,m) 

X6rb 

12 

-0.84570 (0.00112 ) 

1 X 10® 

100 

X6fc 

12 

-5.96339 (0.00157 ) 

1 X 10® 

190 


6 

0.88153 (0.00046 ) 

1 X 10® 

220 

X6fa 

24 

-2.06921 (0.00549 ) 

1 X 10® 

100 

XQjI 

12 

-3.75200 (0.00232 ) 

1 X 10® 

100 


12 

1.64453 (0.00162 ) 

1 X 10® 

200 

XQjM 

6 

3.80656 (0.00138 ) 

1 X 10® 

130 

X6fa 

24 

-0.22601 (0.00143 ) 

1 X 10® 

180 

X<Q3l 

12 

-0.69698 (0.00065 ) 

1 X 10® 

100 

X6fc 

12 

-0.02753 (0.00053 ) 

1 X 10® 

170 


6 

0.72170 (0.00037 ) 

1 X 10® 

100 


46 



TABLE XV: Numerical evaluation of diagrams of Set VI(k) contributing to the muon g-2. The 
superscript e (m) indicates that the lepton loop Ag is the electron (muon) loop. The latter is 
mass-independent so that it is identical with the contribution to the electron g-2. 


Integral 

np 

Value (Error) Sampling per 

No. 

of 



including np 

iteration 

iterations 


10 

50.35921 (0.01998 ) 1 x 10®, 1 x 

10^ 

500 

, 150 

X6kl;^^ 

10 

1.72669 (0.01786 ) 

1 X 

10® 


400 

XQk^c^ 

20 

7.81330 (0.02038 ) 

1 X 

10® 


300 

XQkf 

20 

20.67840 (0.03758 ) 

1 X 

10® 


100 

X6ki^^ 

10 

-0.19466 (0.01045 ) 

1 X 

10® 


300 

XQkf 

10 

1.75890 (0.02374 ) 

1 X 

10® 


230 

XQk^g'^ 

20 

-0.02607 (0.01797 ) 

1 X 

10® 


200 


10 

-0.69054 (0.00750 ) 

1 X 

10® 


100 


10 

15.69736 (0.01595 ) 

1 X 

10® 


200 

XQk^r^ 

10 

-0.56022 (0.00301 ) 

1 X 

10® 


100 

XQk^^^ 

10 

0.30282 (0.00085 ) 

1 X 

10® 


100 

X6kir'^ 

20 

-0.32547 (0.00114 ) 

1 X 

10® 


100 

X6k^P'> 

20 

0.82380 (0.00084 ) 

1 X 

10® 


100 

XQk^r'’ 

10 

-0.17188 (0.00053 ) 

1 X 

10® 


100 

X6k^P'> 

10 

0.30329 (0.00088 ) 

1 X 

10® 


100 

X6k^^'> 

20 

-0.94843 (0.00067 ) 

1 X 

10® 


100 

XQk^^'> 

10 

-0.13877 (0.00018 ) 

1 X 

10® 


100 

XQkf"'> 

10 

1.39510 (0.00069 ) 

1 X 

10® 


100 
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